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Abstract 



Antenna functions describe the infrared singular behaviour of colour-ordered QCD matrix elements 
due to the emission of unresolved partons inside an antenna formed by two hard partons. In this paper, we 
show that antenna functions for hard gluon-gluon pairs can be systematically derived from the effective 
Lagrangian describing Higgs boson decay into gluons, and compute the the infrared structure of the 
colour-ordered Higgs boson decay matrix elements at NLO and NNLO. 



1 Introduction 



Perturbative QCD corrections to exclusive jet observables are at present restricted to the next-to-leading 
order in perturbation theory, which is often insufficient to match the experimental precision on jet production 
reactions [1]. In the recent past, much progress was made to extend these calculations to the next-to-next- 
to-leading order (NNLO) in perturbation theory [2-14] and first results for exclusive NNLO cross sections 
became available recently [15-17]. 

The calculation of NNLO corrections to jet observables requires a method for the extraction of real 
radiation singularities arising from the emission of up to two unresolved (soft or collincar) partons in the 
final state. Several methods have been proposed recently [11] to accomplish the task of constructing so-called 
NNLO subtraction terms. 

In [14], we described the derivation of NNLO subtraction terms for e + e~ — > 2j based on full four- 
parton tree-level and three-parton one-loop matrix elements, which can be integrated analytically over the 
appropriate phase spaces [13]. These NNLO subtraction terms were used subsequently [17] in the calculation 
of the NNLO corrections to one of the colour factors contributing to e+e~ — > 3j. 

Subtraction terms derived from full matrix elements can be viewed as antenna functions, encapsulating 
all singular limits due to unresolved partonic emission between two colour-connected hard partons [6, 18]. 
In particular, process-independent antenna functions describing arbitrary QCD multiparticle processes can 
be directly related to three-parton matrix elements at NLO (one unresolved parton radiating between two 
colour-connected hard partons) and four-parton matrix elements at NNLO (two unresolved partons radiating 
between two colour-connected hard partons). 

QCD calculations of jet observables require three different types of antenna functions, corresponding 
to the different pairs of hard partons forming the antenna: quark-antiquark, quark-gluon and gluon-gluon 
antenna functions. The quark-antiquark antenna functions can be obtained from the e + e~ — ► 2j real 
radiation corrections at NLO and NNLO [14]. In [19], we described how the quark-gluon antenna functions 
could be derived from the purely QCD (i.e. non-supersymmetric) NLO and NNLO corrections to the decay 
of a heavy ncutralino into a masslcss gluino plus partons. It is the purpose of this letter to complete the 
derivation of NLO and NNLO antenna functions by considering the corrections to the decay of a Higgs boson 
into gluons as template for the gluon-gluon antenna functions. 

The Higgs boson coupling to gluons is mediated through massive quark loops, which decouple for large 
quark masses, thus yielding an effective theory containing the interaction of the Higgs field with the gluonic 
field strength tensor [20]. In this effective theory, the Higgs boson decay rate [21] and inclusive Higgs boson 
production cross sections [22-24] were computed to NNLO. Most recently, NNLO results for the exclusive 
Higgs boson production cross section [16] were obtained as well. 

In the following, we will show which individual real radiation processes contribute to the Higgs boson 
decay in the effective theory at NLO and NNLO, and that the real radiation singularities arising at these 
orders precisely match the infrared singularity structure obtained from an infrared factorization formula [25] , 
such that these Higgs boson decay matrix elements can be used to derive the gluon-gluon antenna functions 
at NNLO. 

2 Effective Lagrangian and Feynman rules 

At tree level, the Higgs boson does not couple either to the gluon or to massless quarks. In higher orders 
in perturbation theory, heavy quark loops introduce a coupling between the Higgs boson and gluons. In the 
limit of infinitely massive quarks, these loops give rise to an effective Lagrangian [20] mediating the coupling 
between the scalar Higgs field and the gluon field strength tensor: 

Ant - —HF^Fa,^ . (2.1) 

The coupling A has inverse mass dimension. It can be computed by matching [21,26] the effective theory to 
the full standard model cross sections [27]. 
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The Feynman rules following from this Lagrangian are: 
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The momenta are always incoming. 

In the present context, the value of A is irrelevant, but we do have to take into account that A is 
renormalized in the effective theory. The renormalization constant of A was computed to all orders in [28]; 
it reads 



Zx = 



1 



with 



0o = 



l-0(a,)/e 
11N-2N F 



2?r e 

01 = 





) 2 






V27T, 




e 2 


e 



+ 0{al) 
34iV 3 - 13iV 2 iV F + 3N F 



6 12N 

3 Colour-ordered amplitudes in Higgs boson decay 

The basic process for the decay of a Higgs boson into partons is H(q) — > g(pi)g(p2)- Its amplitude reads 



(2.5) 
(2.6) 



M 



9192 



(3.1) 



The amplitude contains two colour connected (hard) partons which form two antennae, since unresolved par- 
ton emission can take place on both fundamental colour lines connecting the gluons p\ and P 2, as illustrated 
in Figure m 

The squared matrix element, averaged over identical gluons in the final state is 



2 1 9192 
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(3.2) 



Tg g (q 2 ) serves as normalization for antenna functions obtained from higher order corrections to this matrix 
element. 

To demonstrate the cancellation of infrared divergences at NLO, we compute the renormalized one-loop 
QCD correction to the H (q) — > g(pi)g(p2) decay, 
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Figure 1: Colour flow contained in tree level decay H — > gg. Double (single) lines denote adjoint (funda- 
mental) colour indices. 



The infrared poles of this one-loop correction can be expressed in terms of the infrared singularity operator 
[25] 



gg\~>i > 2r(l-e) 
as 
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Voles (T> 2 )) = (g) 4RcI«(e,g 2 ) T> 2 ) . (3.5) 

This expression has to be compared to the ZRel9g (e, q 2 ), which is obtained in the decay of a virtual photon 
into a quark-antiquark pair 7* — > qq at one loop [14] and the factor 4ReJ^(e, q 2 ), which is obtained in the 
decay of a neutralino into a gluino-gluon pair x 99 a t one loop [19]. The factor 4 in (|3.5J) appears since 
the leading order process H — > gg contains two distinct gluon-gluon antennae, just as x ~~ > 99 contains two 
quark-gluon antennae, but in contrast to the single quark-antiquark antenna in 7* — > qq. 

4 NLO antenna functions 

Two different emissions off a gluon-gluon pair appear at NLO: either the emission of an additional gluon or 
the splitting of one gluon into a quark-antiquark pair. In the context of Higgs boson decay, these correspond 
to the tree level processes H — > ggg and H — > gqq. 

The tree level amplitude for H(q) — > g{ P i)g{ P 2)g{Pd,) contains only a single colour structure, f a ^ a 2 a 3 : 

m; iM3 = i\gf aia2a3 M ggg ( Pl ,p 2 , P3 ) . (4.1) 

Squaring the matrix element and dividing by a symmetry factor to account for identical gluons in the final 
state yields 



^ \M° gia2g3 \ 2 = AV (^V 2 - 1) - \M° ggg ( PuP2 , P3 )\ 2 , (4.2) 



with 
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The factor 1/3 in the above equation reflects the fact that the H — > ggg matrix element contains three 
different antenna configurations (corresponding to the three different possibilities of identifying the two hard 
gluons and the one unresolved gluon). The effect of the symmetrization over the three gluons is that these 
three antenna configurations are averaged over. To illustrate the antenna factorization, the leading order 
matrix element (without the symmetrization factor for two identical gluons) is factored out. 

The behaviour of this matrix element in the kinematical limits where one parton becomes unresolved is 
as follows: 
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1. Collinear limits: 
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with z being the momentum fraction of one of the collinear partons and the splitting function 
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(4.4) 
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2. Soft limits: 
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Besides the symmetry factor 1/3 accounting for the average over the three different antenna configurations, 
we observe an overall factor 2, corresponding to the presence of two distinct antenna functions in the basic 
two-parton matrix element, Figure ^ 

To obtain antenna functions describing the emission of an unresolved gluon j off an antenna containing 
two hard gluons i, k, the matrix element l|4.3[) has to be split into three individual antenna configurations. 
Each individual antenna configuration contains only one soft limit. Each collinear g — > gg is split between 
the two antenna configurations appropriate to the two final state gluons involved in the splitting, as discussed 
in [6,18]. 

Integration over the dipole phase space [13] yields 
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The tree level amplitude for H(g) — > g(pi)q{p3)q{pi) contains only a single colour structure T" 1 ^: 



^919394 = iX 9 T rX M gqq(Pl,P3,P4) > 



yielding 



with 



\M° qq J P l,P3,P4)\ 2 



,7V 2 -1 



2 _ 1,2 2tZ -"lA/fO r„ „ „.M2 

919394 I n i> r 



M° (pi,p 3 ,p 4 )| 2 , 



(1-0 2 



(fl3 + Sl4) 
S34 



- 2 



S34 



(4.6) 
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The only singular configuration contained in this matrix element is the collinear quark- antiquark limit, which 
is as follows: 

' ?4 2 (47ra s ) T g ° g ( SlM )— P a ^ qq (z) , 
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with the collinear splitting function 
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The factor 2 arises from the fact that two gluon-gluon antennae are contained in the matrix element l|3.1|l 
as above. 
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Integration over the dipole phase space [13] and summing over final state quark flavours yields 

T g \ q (q 2 ) ee [d* D , gq ,J2\Mi q3Cl f 
J q 

= (g)^W(s 2 r 
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Summing over both three parton final states, we find 

™e« (T s y<z 2 ) + T s y, 2 )) = - (g) 4Re7« (e , , 2 ) T> 2 ) 
such that the NLO corrected Higgs boson decay rate into partons is finite: 

Voles (T/ ff ( g 2 )) + Voles {^(q 2 ) + T g \ q (q 2 )) = . 
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We recover the finite NLO contribution to the Higgs boson decay rate into partons in the effective theory as 
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which is in agreement with [21]. 



5 Structure of NNLO antenna functions 

In the NNLO calculation of jet observables, two different types of antenna functions are required: (a) the 
one-loop correction to the three-parton antenna functions which appeared at NLO in tree-level form, and 
(b) the tree-level four-parton antenna functions. In this section, we present all Higgs boson decay matrix 
elements needed for the derivation of these antenna functions, and demonstrate that these matrix elements 
contain the same infrared singularities as processes involving final state emission off a gluon-gluon antenna. 

The renormalized one-loop corrections to the three-parton antenna functions have the same colour struc- 
ture as their tree level counterparts listed above. To expose the infrared structure of the resulting one-loop 
matrix elements, they are integrated over the corresponding dipole phase space [13], yielding 
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Figure 2: Colour flow contained in the colour ordered amplitude M g (pi,pi,pj,pk) contributing to the tree 
level decay H — > gggg- 
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Figure 3: Colour flow contained in the colour ordered amplitude Mg gqq (pi, P2, P3, Pa) contributing to the tree 
level decay H — > ggqq. 
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Three different four-parton final states appear in the gluon-gluon antenna functions at NNLO: gggg, qqgg 
and qqq'q' . In contrast to the tree level three-parton Higgs boson decay matrix elements, which contained only 
one non-trivial colour ordering each, these four-parton matrix elements all contain several colour-orderings. 

The amplitude for H(q) — > g(pi)g{p2)g(P3)g{PA) can then be expressed as sum over the permutations of 
the gluon colour indices: 
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(5.3) 



where the sum runs over all six permutations of three of the gluon colour indices, thus excluding any 
configurations which can be related by cyclic permutations of all four colour indices. Its colour flow is 
illustrated in Figure |2 

The resulting squared matrix element, averaged over identical final state gluon permutations is 
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It should be noted that this squared matrix element contains only the leading colour term obtained from 
the squares of the individual colour-ordered amplitudes, as expected in the colour ordered formulation for a 
process with four gluons [29,30]. 

The tree level amplitude for H(q) — > g{pi)g{p2)q{p3)q(pi) contains two colour structures, 
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The squared matrix element, averaged over identical gluons in the final state and summed over quark 
flavours, reads 
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Finally, the tree level amplitude for H(q) — > q(pi)q{p2)q' '{pz)<l '{Pi) contains only a single colour structure, 
but can contain two flavour structures in the case of identical quark flavours q = q'\ 
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The squared matrix element, summed over quark flavours and averaged over identical configurations 
becomes 
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The four-parton tree-level Higgs boson decay matrix elements can be integrated over the tripole phase 
space [13], thus making their infrared singularity structure explicit, 
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The sum of all NNLO subtraction terms yields the following infrared pole structure, which can be 
expressed in terms of NNLO infrared singularity operators [25], 
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where /3o is the first term of the QCD /3-function (|2.6|) and the constant K 
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The above structure is to be compared with the renormalized purely virtual NNLO corrections (two-loop 
times tree plus one-loop self- interference) , which were first computed by Harlander [31]: 
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It can be seen that the poles of the real radiation terms (|5.12|l cancel the poles of the purely virtual 
corrections: 
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The infrared singularity structure Ij5.13|) corresponds to the NNLO corrections to a tree level process con- 
taining two gluon-gluon antenna functions, as is the case for the Higgs boson decay. 
An important check on our results is that the sum of all NNLO contributions 
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agrees with the NNLO correction to the total Higgs boson decay rate into partons in the effective theory [21]. 

Equations l|5.13[l and (|5.16|) demonstrate that the NNLO three and four parton contributions to Higgs 
boson decay into massless partons display the same singularity structure as final state observables containing 
adjacent gluon-gluon pairs. It is therefore possible to derive colour-ordered gluon-gluon antenna functions 
from the Higgs boson decay matrix elements obtained here using the effective Lagrangian density (|2.1|) . 



6 Conclusions and Outlook 

QCD antenna functions describe the behaviour of QCD matrix elements in their infrared singular limits, 
corresponding to soft or collinear parton emission. They are constructed so that they describe all singular 
limits arising from emission of unresolved partons in between the two colour-connected hard partons that 
define the antenna. The quark- antiquark antenna functions are directly related to the physical matrix 
elements for 7* — » gg+partons. We demonstrated in a previous paper [19] that quark-gluon antenna functions 
could be obtained from an effective Lagrangian density describing ncutralino decay into a gluino and other 
partons. Besides quark-antiquark and quark-gluon antenna functions, QCD calculations also require gluon- 
gluon antenna functions. In this paper, we showed that gluon-gluon antenna functions can be obtained from 
physical Higgs boson decay matrix elements into partons, arising in an effective theory coupling the Higgs 
field to the gluonic field strength tensor. 

We demonstrated that the physical Higgs boson decay matrix elements reproduce the singular structure 
of QCD gluon-gluon antenna functions at NLO and NNLO. We extracted the infrared structure for decay 
kinematics, as required for jet observables without partons in the initial state. By analytic continuation, the 
matrix elements derived here can also be continued to production (leading order process contains partons 
only in the initial state) or scattering (leading order process contains partons in initial and final state) 
kinematics, where they have to be integrated over the appropriate phase spaces. The phase space integrals 
for production kinematics were derived in [24] , such that the antenna subtraction terms for this kinematical 
situation can in principle be derived. 

With this and two preceeding papers [14,19], we demonstrated that all QCD antenna functions can 
be derived (as opposed to constructed) from physical matrix elements: quark-antiquark antennae from the 
decay of a virtual photon into partons, quark-gluon antennae from neutralino decay into gluino plus partons 
and finally gluon-gluon antennae from Higgs boson decay into partons. The NNLO antenna subtraction 
functions obtained through this procedure will be reported in a subsequent publication [32]. 
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